An Attempt to Probe the Radio Jet Collimation Regions in 
NGC4278, NGC4374 (M84), and NGC 6166 
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ABSTRACT 

NRAO Very Long Baseline Array (VLBA) observations of NGC 4278, NGC 4374 (M84), 
NGC 6166, and M87 (NGC 4486) have been made at 43 GHz in an effort to image the jet coUi- 
mation region. This is the first attempt to image the first three sources at 43 GHz using Very 
Long Baseline Interferometry (VLBI) techniques. These three sources were chosen because their 
estimated black hole mass and distance implied a Schwarzschild radius with large angular size, 
giving hope that the jet collimation regions could be studied. Phase referencing was utilized for 
the three sources because of their expected low flux densities. M87 was chosen as the calibrator 
for NGC 4374 because it satisfied the phase referencing requirements: nearby to the source and 
sufficiently strong. Having observed M87 for a long integration time, we have detected its sub- 
parsec jet, allowing us to confirm previous high resolution observations made by Junor, Biretta 
& Livio, who have indicated that a wide opening angle was seen near the base of the jet. Phase 
referencing successfully improved our image sensitivity, yielding detections and providing accu- 
rate positions for NGC 4278, NGC 4374 and NGC 6166. These sources are point dominated, but 
show suggestions of extended structure in the direction of the large-scale jets. However, higher 
sensitivity will be required to study their sub-parsec jet structure. 

Subject headings: galaxies: active - galaxies: elliptical and lenticular, cD - galaxies: individual 
(NGC 4278, NGC 4374, NGC 4486, NGC 6166) - galaxies: jets - radio continuum: galaxies 



1. INTRODUCTION 

The high resolution capabilities of Very Long 
Baseline Interferometry (VLBI) have given as- 
tronomers the ability to study jets in the inner 
parsec of active galactic nuclei (AGN), leading to 
the possibility of testing jet collimation theories 
and magnetohydrodynamic (MHD) simulations. 
MHD simulations have been able to explain the 
formation of AGN jets by magnetic fields thread- 
ing through an accretion disk around a super- 
massive black hole. The fields thread out from 
the disk in a rotating helical coil, accelerating and 
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guiding the plasma away from the disk. The re- 
sulting jet structure is expected to show a large 
opening angle near the base and tight collimation 
on larger scales (Meier, Koide, & Uchida 2001). 
An example where the large opening angle at the 
base of the jet may have been observed is the ellip- 
tical galaxy, M87 (Junor, Biretta, & Livio 1999). 
They reported that the jet must form within 30 
Schwarzschild radii (i?s) of the black hole and 
the collimation region extends to beyond 100 Rs 
(Biretta, Junor, & Livio 2002). M87's close prox- 
imity, large black hole mass, and bright jet make 
it arguably the best candidate to study the jet col- 
limation region. 

Motivated by M87, we began an effort to ob- 
serve other sources in which the jet collimation re- 
gion might be probed. We selected sources based 
on the expected angular size of the Schwarzschild 
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radius [Rs = {2GM,)/c^], since that is likely to 
set the scale of the collimation region. Hence, 
nearby galaxies with large black hole mass were 
favored. Wc further selected sources based on 
the presence of detectable radio emission. Using 
black hole masses available in mid 2000, the most 
promising sources were NGC 4278, NGC 4374 and 
NGC 6166. The estimated black hole masses avail- 
able then were 15, 18 and 290x10^ Mq, respec- 
tively (Gebhardt et al. 2000; Magorrian et al. 
1998). NGC 6166 is relatively distant, but what 
turned out to be an excessively large black hole 
mass estimate qualified it as part of our sample. 

A recent study finds that the correlation be- 
tween the mass of the central black hole in a galaxy 
and the velocity dispersion (M, — a relation) is 

log (M./Mq) = a + /Jlog (a/ao) (1) 

where a = 8.13 ± 0.06, /3 = 4.02 ± 0.32, M. is the 

black hole mass, a is the central velocity dispersion 
in units of km s~^ and <Jo is a reference velocity 
taken to be 200 km s"^ (Tremaine et al. 2002). 
For our three target sources, the black hole masses 
computed from Equation 1 are substantially lower 
than those used for the selection process. 

Current distance and black hole mass estimates 
are provided in Table 1 for our target sources and 
M87. The first and second columns list the NGC 
and 3C catalog numbers of the sources; Messier 
catalog names are in parentheses. Column 3 pro- 
vides the computed distances D in units of Mpc to 
these sources using surface brightness fluctuation 
(SBF) distance moduli, and Column 4 gives the 
references for these values. Black hole masses are 
tabulated (Column 5) in units of 10^ Mq. These 
numbers were computed using Equation 1 and the 
velocity dispersions from the references (Column 
6). For M87; the black hole mass is as reported 
by Tremaine et al. (2002). Using the values in 
Columns 3 and 5, the inferred angular sizes in 
micro- arcseconds (/ias) of the Schwarzschild radius 
is provided in Column 7. 

Section 2 will describe our observing tech- 
niques, the calibration process and our imaging 
approach. Section 3 will briefly introduce charac- 
teristics of each source, including previous related 
observations, and discuss the implications inferred 
from our new VLBA images. Finally, Section 4 
will summarize our results, and remark on the 
future prospects for probing the jet collimation 



regions. 

2. OBSERVATIONS, CALIBRATION, 
AND IMAGING 

Observations with the 10 VLBA (Napier et 
al. 1994) antennas of NGC 4374 (also M87), 
NGC 6166 and NGC 4278 were made on 2001 Oc- 
tober 12, 2002 February 14, and 2002 Septem- 
ber 19, respectively. These observations were 
made at 43.127 GHz using 32 MHz bandwidth 
in each polarization and 2 bit sampling. The res- 
olution was typically 0.2 by 0.4 milli-arcseconds 
(mas). Since the target sources were faint, phase- 
referencing techniques were used to improve im- 
age sensitivity. Suitable phase calibrators (peak 
flux density greater than 100 mJy beam~^ and 
within 2° from target sources) were J1221-I-2813 
(W Comae) for NGC 4278, M87 for NGC 4374, 
and J1635-I-3808 for NGC 6166. For the first ob- 
servation, NGC 4374 and M87 were each observed 
with short, scheduled scans of 15 seconds, using 
the "nodding" technique. Analysis of that obser- 
vation showed that only a few seconds of on-source 
time were obtained for each scan, and the coher- 
ence time was adequate to allow a longer switch- 
ing time. Therefore, the other two observations 
were made with scans of 30 seconds on target and 
20 seconds on calibrator, for a total switching time 
of 50 seconds. The approximate total integration 
times were 200 minutes NGC 4278, 91 minutes for 
NGC 4374, 85 minutes for M87, and 141 minutes 
for NGC 6166. 

In addition to our target and calibration 
sources, we also observed fringe finders and check 
sources. The fringe finders, 3C345 and 3C279, 
were used to make delay corrections. Check 
sources were included to examine the quality of 
the phase referencing. The check sources were 
J1218-hll05, J1217+3007 and 3C345 while ob- 
serving NGC 4278, NGC 4374 and NGC 6166, re- 
spectively. 

We followed the standard VLBA data reduc- 
tion procedure using the NRAO Astronomical Im- 
age Processing System (AIPS) as described in Ap- 
pendix C of the AIPS Cookbook. Delays were cor- 
rected using fringe fitting with the fringe finders. 
They were sufficiently constant that no further 
fringe fitting was done. Phase calibrators were it- 
eratively imaged and self-calibrated. The images 
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were then used as models to determine the am- 
plitude and phase gains for the antennas. These 
gains were applied to the target and check sources 
(phase referencing) and images were made. The 
target sources proved to be strong enough to al- 
low slowly varying gain errors to be corrected and 
improved images to be made using self-calibration 
with long (10, 30 or 60 minutes) solution inter- 
vals. Saint Croix was excluded for the phase refer- 
enced sources because of rapid phase fluctuations, 
but was used for the self-calibrated images of the 
sources observed on 2001 October 12. 

Observational results on all sources, excluding 
fringe finders, are provided in Table 2. We have 
included both phase-referenced and self-calibrated 
flux densities to indicate the quality of the phase 
referencing. Columns 1 and 2 provide the name 
of the source and what it was used for during the 
observation. 'Calib', 'Target' or 'Check' indicates 
that the source was used either as a calibrator, a 
target or a check source. Columns 3 and 4 list the 
right ascension and its one-sigma error with the 
declination and its one-sigma error below it. The 
coordinates and errors for the calibrators were ob- 
tained from the International Celestial Reference 
Frame Extension 1 (ICRF-Ext.l; lERS 1999; Ma 
et al. 1998); the positions for the other sources 
were measured from our images. Errors account 
for the error of the calibrator position (the domi- 
nate source of error), an assumed 0.1 mas poten- 
tial error due to differential geometry and atmo- 
spheric model errors between calibrator and target 
or check source, and errors in the measurement of 
position from the images. Column 5 reports the 
angular offset (9) in degrees from the calibrator. 
Columns 8 through 13 are the peak flux density 
per beam (Sp) of the source in mJy beam^^, the 
integrated flux density {Si) in mJy, and the off- 
source image noise (RMS) in mJy beam~^ for 
both the phase-referenced and the self-calibrated 
images. To obtain the integrated flux for our tar- 
get sources and M87, we used a region that is 
2.5 mas x 1 mas for NGC4278, 2.8 mas x 1 mas 
for NGC4374, 2.5 mas x 5 mas for M87, and 
1.5 mas X 1.5 mas for NGC6166. Phase- 
referenced and self-calibrated images have been 
convolved to the same beam size and position 
angle listed in Column 5 and 6. For 3C345, ob- 
served on 2002 February 14, only the portions of 
the data scheduled as a check source were included 



in the phase referencing. The self-calibrated val- 
ues include all data. The uncertainty of the self- 
calibrated flux densities is difScult to determine, 
but wc estimate it to be about 10%. It depends on 
the quality of the VLBA a priori gain calibration, 
the ability of the phase referencing to remove short 
term fluctuations for weak sources, the quality of 
the self-calibration, and the image noise. In all 
cases, the flux scale (gain normalization) was set 
by observations above 30° elevation at a subset 
of antennas that were well-behaved and had good 
weather on the day of the observation. 

The 100-y^as potential error in the positions due 
to differential model errors is approximate and 
was derived in two ways. First, the raw phases 
showed approximately 10 turns of phase across the 
day, which is equivalent to across the sky on the 
longest baselines. Scaling this to an angular er- 
ror, and decreasing it by the ratio of the calibra- 
tor/source separation to a radian, gives an esti- 
mate of about 0.1 mas. Second, Pradel, Char- 
lot, & Lestrade (2003) have done a study of the 
errors in differential VLBI astrometry, suggesting 
errors of approximately 0.1 mas for our observa- 
tions (Lestrade, private communication). Two of 
our check sources (J1218-hll05, 3C345) and two 
of our target sources (NGC4278, NGC4374) are 
included in the ICRF-Ext.l. In all cases, the 
calibrator-target or calibrator-check source differ- 
ential position measured from our images agrees 
with the differential positions from the ICRF- 
Ext.l to better than one- third to one- half of the 
reported ICRF-Ext.l errors. This confirms that 
our measured positions are reasonable and implies 
that both our differential positions, and the rela- 
tive positions of close pairs of sources in the ICRF- 
Ext.l, are significantly more accurate than the ab- 
solute position errors in the ICRF-Ext.l. 

3. DISCUSSION 

3.1. NGC 4278, NGC 4374 and NGC 6166 

Our VLBA images of NGC 4278 (Figure la), 
NGC 4374 (Figure lb), and NGC 6166 (Figure Id) 
are currently the highest resolution images avail- 
able for these galaxies. Continuum emission has 
been detected from the cores at 43 GHz with 
FWHM resolutions of about 400 /xas x 200 jias. 
While these VLBA images do not clearly reveal 
any jet structures, they do provide indications of 



3 



low-lcvcl emission that may be associated with 
jets traced at lower frequencies and lower resolu- 
tions. Before discussing each galaxy below, we 
note that the photometry of the cores at sub-mas 
resolution (Table 2) may help constrain spectral 
energy distributions of any advection-dominated 
accretion flows in these galaxies (Di Mattco, Car- 
illi, & Fabian 2001; Di Matteo, Johnstone, Allen, 
& Fabian 2001; Verdoes Kleijn, Baum, de Zeeuw, 
& O'Dea 2002). 

3.1.1. NGC4278 

Continuum studies of this elliptical galaxy at 
arcsecond resolution show dominance by a flat- 
spectrum, time-variable component at frequencies 
of 1-15 GHz (Wrobel & Heeschen 1984, 1991; Na- 
gar, Wilson, & Falcke 2001; Nagar, Falcke, Wil- 
son, & Ulvestad 2002). At a FWHM resolution of 
3" (230 pc), this component is linearly polarized 
at 8.4 GHz by 0.34 ± 0.06 percent with an elec- 
tric vector PA of 16 ± 7° (Bower, Falcke, & Mel- 
lon 2002). Assuming Faraday rotation effects are 
minimal at this frequency, this polarized emission 
probably arises in the diffuse structure, mainly 
to the south and southeast, traced in images at 
8.4 GHz at a FWHM resolution of 200 mas (16 pc; 
Wilkinson et al. 1998) and at 5 GHz at FWHM 
resolutions of 400 mas (31 pc; Wrobel & Heeschen 
1984) and 2.5 mas (0.2 pc; Falcke, Nagar, Wil- 
son, & Ulvestad 2000; Giovannini et al. 2001). 
In Figure la, the arrow indicates the direction to 
the jet component offset from the 5 GHz core by 
10 mas (0.8 pc), as measured with the VLBA by 
Falcke, Nagar, Wilson, & Ulvestad (2000). While 
our VLBA image of NGC 4278 at 43 GHz shows 
hints of low-level emission mainly to the north, 
that emission is too close to the noise to be con- 
sidered reliable. Finally, we note that our VLBA 
image at 43 GHz recovers less than one half of 
the 43-GHz flux density of the unresolved VLA 
core, in 2000 February, at a FWHM resolution 
of about 500 mas (39 pc; Di Matteo, Carilli, & 
Fabian 2001). Such behavior is expected because 
the VLBA is insensitive to mas-scale emission at 
43 GHz, known to be present at lower frequen- 
cies. Still, some of this difference could originate 
through variability of the sub-mas core between 
2000 February and 2002 September. 



3.1.2. NGC 4374 (M84) 

This elliptical galaxy hosts a well-studied 
Fanaroff-Riley class I (FR-I) radio continuum 
source (Laing fc Bridle 1987). Images at 5 GHz 
at a FWHM resolution of 500 mas (44 pc) show 
two large-scale jets that are initially asymmetric, 
with the northern main jet being brighter than 
the southern counter] et. Images at 1.7 GHz at 
FWHM resolutions of 11 mas (1 pc; Xu et al. 

2000) and of 6 mas (0.5 pc; Giovannini et al. 

2001) show extensions to the north. A VLBA im- 
age at 5 GHz at a FWHM resolution of 2.7 mas 
(0.2 pc) shows a core plus two components to the 
north and one to the south, offset in the direc- 
tions indicated by the arrows in Figure lb (Nagar, 
Falcke, Wilson, & Ulvestad 2002). The available 
ICRF-Ext.l position and 1-a error for NGC 4374 
are marked by a cross in Figure lb (lERS 1999). 
Our VLBA image of NGC 4374 at 43 GHz shows 
that the core appears to be slightly extended to 
the north. Unfortunately, the image sensitivity 
for this galaxy was reduced compared to the other 
two target sources, due to our less efficient ob- 
serving strategy. For this reason, and because it 
is strong enough not to require phase referencing, 
NGC 4374 is a prime choice for future observations 
with significantly improved sensitivity. 

3.1.3. NGC6166 (3C338) 

This cD galaxy with multiple nuclei hosts an 
FR-I source (Giovannini et al. 1998, plus refer- 
ences therein). VLBI images between 1.7 and 
8.4 GHz with FWHM resolutions as high as 
1.2 mas (0.7 pc) reveal two-sided parsec struc- 
ture in the direction of the arrows in Figure Id. 
Our VLBA image of NGC 6166 at 43 GHz shows 
slight elongations along those directions. It is un- 
fortunate that the black-hole mass estimate has 
been revised downward by a factor of 40, offering 
no real prospect for probing the jet coUimation 
region with this VLBA observation. 

3.2. M87 

This elliptical galaxy hosts an FR-I source 
whose main jet is especially well-studied on arc- 
second and milli- arcsecond scales (Biretta, Junor, 
& Livio 2002, plus references therein). Global 
VLBI observations at 43 GHz offer sub-mas reso- 
lution and indicate that the initial jet coUimation 
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region extends to at least 0.5 mas (Junor, Biretta, 
& Livio 1999), corresponding to 0.04 pc or 140 Ra 
for the black hole mass and galaxy distance in Ta- 
ble 1. Figure Ic shows our VLB A image of M87 
at 43 GHz, with a FWHM resolution of 390 fx&s x 
170 /xas, which independently confirms the Junor 
et al. result of coUimation over the inner 1 mas. 
Moreover, jet emission is detected out to about 
3 mas (0.2 pc or 820 Rs) along the direction of 
the large-scale jet indicated by the arrow in Fig- 
ure Ic. One difference between our image and 
that of Junor et al. is the presence of a weak 
counterjct in our image. Detection of a coimter- 
jet would be very interesting, if real. But such 
structures are also a natural consequence of am- 
plitude calibration errors. We attempted, through 
manipulation of the region in which the CLEAN 
algorithm was allowed to place components during 
the sclf-calibration/imaging iterations, to make 
images without the counterjet. While success- 
ful, the counterjet tended to reappear when the 
CLEAN boxes were opened up to allow it. In the 
end, despite considerable testing, the counterjet 
detection can only be considered suggestive, not 
reliable. 

3.3. Other sources 

One of our calibrators, J1221-I-2813 (W Co- 
mae), showed extended structure with a compo- 
nent 3 mas southeast from the bright core. The 
structure is consistent with that seen at 22 GHz 
by Wiik, Valtaoja, & Leppanen (2001). Our 
observation of 3C345 shows a knot 1 mas west 
from the radio core with additional extended emis- 
sion. This is consistent with 22 GHz observa- 
tions made by Ros, Zensus, & Lobanov (2000). 
J1218-I-1105 appears to be a compact source at 
43 GHz. The check source for our third obser- 
vation (J1217-I-3007) is seen mostly as a compact 
source with low- level emission to the west. Our 
brightest calibrator, J1635+3808 has a minor jet 
to the west, appearing "thumb- like" . 

4. CONCLUSION 

VLBI observations at 43 GHz have shown that 
M87 jet has a wide opening angle in the inner 
100 Rg, providing important clues to jet coUima- 
tion mechanisms (Biretta, Junor, & Livio 2002). 
The VLBA observations reported here were the 



first from an effort to identify other sources in 
which the structure of the base of the jet might 
be observed. The target sources were chosen to 
have radio continuum emission and large, angular- 
size black holes. Each of the three target sources, 
NGC4278, NGC4374, and NGC6166 (3C338), 
had point-dominated structure with only hints of 
the beginning of the jet. The sensitivity of the 
observations proved to be inadequate to deter- 
mine the structure of the jet bases. M87 was 
also observed as the calibrator for NGC 4374; the 
resulting image agrees with the results of Junor, 
Biretta, & Livio (1999). Future observations will 
concentrate on improving sensitivity and finding 
new sources. 

The observations reported here involved phase 
referencing at 43 GHz with the VLBA. The tech- 
nique worked well with switching angles between 
calibrator and target of 1-2° and switching time 
intervals of 50 seconds (20 seconds on calibrator, 
30 seconds on target). A 30 second switching time 
(i.e. 15 seconds per scan) was tried, but gave 
unreasonably short on-source integration and was 
not required by the coherence time. The methods 
used here allow detection of sources as weak as 
2 mJy at 43 GHz with ~ 8-hour observations. 

This research has made use of the NASA/IPAC 
Extragalactic Database, which is operated by the 
Jet Propulsion Laboratory. Caltcch, under con- 
tract with the National Aeronautics and Space 
Administration, and NASA Astrophysics Data 
System. C. Ly acknowledges support at the Na- 
tional Radio Astronomy Observatory by the Re- 
search Experiences for Undergraduates program 
of the National Science Foundation. 
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Table 1 

Distance and Black Hole Mass Estial4tes 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 






D 




M. 






NGC 


3C 


(Mpc) 


Rcf. 


(10«Mq) 


Rcf. 


(/ias) 


4278 • • • 




16.1 


1 


3.4lj 


2 


0.42 


4374(M84) •• 


■ 272.1 


18.4 


1 


6.311:? 


2 


0.68 


4486(M87) 


• 274 


16.1 


1 


3011° 


3 


3.67 


6166 ••• 


338 


125 


4 


7.111:1 


2 


0.11 



References. — (1) Tonry et al. 2001; (2) Ravindranath, Ho, & 
Filippenko 2002; (3) Tremaine et al. 2002; (4): Jensen et al. 2001 
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Table 2 

vlba astrometry and photometry at 43 ghz 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 

RA (a) Phasc-rcforonced Self-calibrated 

Source Source DEC (5) Aq (ms) 6 FWHM PA ^ RMS Sl ^ RMS~ 

name usage (J2000) A(5 (mas) (°) (mas x mas) (°) (mjy) (mjy) (mjy) (mJy) (mJy) (mJy) 

Observed on 2001 October 12 

NGC4486 ••• Calib 12 30 49.423381 0.017 ••• .39x.l7 168 ••• ••• ••• 610 1300 0.49 

+12 23 28.04393 0.25 

NGC4374 ••• Target 12 25 03.743336 0.018 1.5 .47x.l9 167 45 80 0.86 71 100 0.56 

+12 53 13.13934 0.27 

J1218+1105 ••• Check 12 18 26.092289 0.018 3.3 .66x.23 167 41 69 3.2 69 95 1.8 
+11 05 05.26293 0.27 

Observed on 2002 February 14 

J1635+3808 ••• Calib 16 35 15.492973 0.055 ••• .38x.21 159 5100 5600 1.6 

+38 08 04.50060 0.66 

NGC6166 ••• Target 16 28 38.244693 0.056 1.9 .40x.22 163 14 22 0.43 15 23 0.44 

+39 33 04.23414 0.67 

3C345 ■■■ Check 16 42 58.809985 0.056 2.2 .41x.21 156 1400 .3900 39 1500 4600 2.8 

+39 hS :',(). 99 1:!2 O.GT 

Observed on 2002 September 19 

J1221+2813 ••• Calib 12 21 31.690515 0.019 ••• .42x.l7 158 ••• ••• ••• 88 280 0.5 

+28 13 58.50026 0.29 

NGC4278 ••• Target 12 20 06.825429 0.021 1.1 .40x.l7 159 6.2 19 0.32 8.2 18 0.32 

+29 16 50.71418 0.31 
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Fig. 1. — VLBA images of Stokes I at 43 GHz. Contour levels are scaled from the lowest level by -2.8, -2, -1, 
1, 2, 2.8, 4, 5.7, 8 and multiples of \/2 thereafter, (la) NGC 4278: lowest level is .5 mJy beam~^; grayscale 
levels range from -1.4 to 8.2 mJy bcam"-'^. (lb) NGC 4374: lowest level is 1 mJy bcam~^; grayscale levels 
range from -1.9 to 20 mJy beam~^. {Icj M87: lowest level is 1 mJy beam~^; grayscale levels range from 
-2.5 to 25 mJy beam~^. (Id) NGC 6166: lowest level is 0.5 mJy beam~^; grayscale levels range from -2.0 
to 15 mJy beam"'^. The arrows indicate the direction of the jets seen in low frequency VLBI observations. 
The position offsets are relative to the a priori positions used in correlation, except for NGC 6166 where 
position corrections were made. The cross in lb marks the ICRF-Ext.l position. The beam sizes and peak 
flux densities (self-calibrated) can be found in Table 2. The apparent counterjet in M87 is still questionable 
since it may be an imaging artifact. 9 



